Quantum phase transition in the Frenkel-Kontorova chain: 
from pinned instanton glass to sliding phonon gas 
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We study analytically and numerically the one-dimensional quantum Frenkel-Kontorova chain 
in the regime when the classical model is located in the pinned phase characterized by the gaped 
phonon excitations and devil's staircase. By extensive quantum Monte Carlo simulations we show 
that for the effective Planck constant h~ smaller than the critical value h c the quantum chain is 
in the pinned instanton glass phase. In this phase the elementary excitations have two branches: 
phonons, separated from zero energy by a finite gap, and instantons which have an exponentially 
small excitation energy. At h = h c the quantum phase transition takes place and for h > h c the 
pinned instanton glass is transformed into the sliding phonon gas with gapless phonon excitations. 
This transition is accompanied by the divergence of the spatial correlation length and appearence 
of sliding modes at h~ > h c . 

PACS numbers: 05.45.-a, 63.70.+h, 61.44.Fw 



I. INTRODUCTION 

The Frenkel-Kontorova (FK) model [Q describes a one- 
dimensional chain of atoms /particles with harmonic cou- 
plings placed in a periodic potential. This model was 
introduced more than sixtyyears ago with the aim to 
study crystal dislocations Jj], It was also success- 
fully applied later to the description of commensurate- 
incommensurate phase transitions epitaxial mono- 
layers on the crystal surface ||, ionic conductors and 
glassy materials Jjj, |(], Q| and, more recently, to charge- 
density waves H and dry friction M, |l(|. Despite the 
fact that the relevant phenomena are at the atomic scale, 
all these works are based essentially on the classical ap- 
proach. The first study of quantum effects was done 
twelve years ago Q, |l2| [l3| with the attempt to under- 
stand the highly nontrivial quantum ground state of the 
Frenkel-Kontorova model in the regime when the clas- 
sical ground state is characterized by the fractal "devil 
staircase" ||. These studies were extended in |fLl| [T^| 
and later the quantum dynamics at different values of 
quantum parameter H was studied in [|l6|, |l7], |l8[ |l9|, po[ . 

The physical properties even of the classical FK model 
are very rich and nontrivial. In 1978 Aubry discovered 
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|| a new type of ground state which has fractal prop- 
erties known as "devil's staircase". In fact, the equilib- 
rium positions of atoms in the FK chain are described 
H by the well known Chirikov standard map J2lj], which 
describes generic properties of chaotic Hamiltonian dy- 
namics. The density of particles in the equilibrium state 
determines the rotation number of the invariant curves of 
the map, while the amplitude of the periodic potential in 
the FK model gives the value of the dimensionless param- 
eter K. According to the known properties of the map, it 
follows that there exist two phases of the chain, the "slid- 
ing" phase and the "pinned" phase. Indeed at K < K c 
the Kolmogorov-Arnold-Moser (KAM) curves are smooth 
and the chain can easily slide along the potential. This 
implies the existence of a zero phonon mode. In con- 
trast, at K > K c the KAM curves are destroyed and re- 
placed by an invariant Cantor set which is called cantorus 
H |2^, |2ij, . In this "pinned" phase the chain cannot 

slide being kept by a finite Peierls-Nabarro barrier and 
the phonon spectrum is separated from zero by a finite 
gap. In this paper we consider only the gaped/pinned 
phase with K > K c . 

It is important to stress that in the pinned phase, be- 
sides the equilibrium state of minimal energy, there exists 
a lot of other equilibrium configurations, corresponding 
to local minima of the potential with energies very close 
to the minimal energy, i.e. the energy of the ground state 
The total amount of these states (configu- 
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rational excitations) grows exponentially both with the 
length L of the chain, as well as with the parameter K 
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pgfl . Moreover, a great number of them is practically de- 
generate since their energy separation from the ground 
state is exponentially small. 

In the classical limit all these conflgurational excita- 
tions are stable, while in the quantum case they become 
metastable due to tunneling between different exponen- 
tially degenerate minima of the potential. As a result 
one may expect that the true eigenstates of the Hamilto- 
nian, including the ground state, are built as superposi- 
tions of classical conflgurational states. Such a nontrivial 
structure of the quantum ground state of the chain has a 
profound analogy with the famous vacuum of Quantum 
Chromodynamics (QCD) j2^] in which tunneling tran- 
sitions between different, practically degenerate, states 
are known as "instantons". This analogy between the 
two problems is very useful and implies that many of 
the methods developed in lattice QCD studies may be 
applicable to quantum FK model. 

Dynamical low energy excitations in the classical FK 
chain consist only of phonon modes describing small vi- 
brations around a classical minimum of the chain po- 
tential energy. In the quasiclassical regime, where the 
effective dimensionless Planck constant H is very small, 
the quantization of the FK chain can be reduced to a 
quantization of phonon modes only. Indeed the time of 
tunneling between different minima of the potential en- 
ergy is exponentially large, compared to periods of the 
vibrations, and the phonon modes are decoupled from 
the tunneling modes (instantons) . This reminds the sit- 
uation in the QCD |2j| where the quasiclassical regime 
for instantons appears at small distances on which in- 
stantons are also decoupled from other excitations like 
quarks and gluons. 

In this regime the tunneling transitions are very slow, 
instantons are local and frozen in space. Due to exponen- 
tial degeneracy of chain configurations the model reveals 
a glass-like structure of instantons randomly distributed 
along the chain. We will call this phase the "instanton 
glass". 

One may expect that the phase structure changes sig- 
nificantly with the increase of h when the tunneling time 
t ~ exp(Sin/K) (Sin is the instanton action) becomes 
comparable to the inverse frequency of phonons. In this 
case phonon and instanton excitations can strongly influ- 
ence each other due to anharmonicity of the periodical 
potential. This may lead to quantum melting of the in- 
stanton glass phase and transition to another phase which 
appears above some critical value % > h c . 

The existence of another quantum phase can be ar- 
gued in the following way. At sufficiently large h > h c 
the kinetic energy of quantum particle Ek ~ ti 2 /2m(Ax) 2 
starts to exceed the height U of the periodic potential in 
which the chain is placed (here Ax is the mean particle 
separation in the chain which is comparable with the pe- 
riod of the potential, m is the particle mass). The condi- 
tion E], ~ U gives h c ~ Axy/ mU . As a result, for h> h c 
the chain turns from the pinned to the sliding phase. In 
this regime tunneling is replaced by direct propagation 



above the barrier. Hence, the instanton-like motion is 
replaced by phonon-like motion corresponding to a new 
phase with gapless phonon excitations. This regime can 
be considered as sliding quantum phase similar to classi- 
cal sliding regime at K < K c p^ , B4| . Qualitative differ- 
ent types of the behavior of quantum FK model at small 
and large K values were already seen in the first numer- 
ical studies l^lj |l2|| . In summary, at h = H c one may 
expect a quantum phase transition, with qualitative re- 
arrangement of the spectrum of elementary quasiparticle 
excitations. 

In this work, devoted to a comprehensive numerical 
study of the transition we have just outlined, we present 
a complete picture of low energy excitations in the quan- 
tum FK chain. We also find that this transition is highly 
nontrivial and the model reveals a non-analytical behav- 
ior which can partially explain the failure of simple ana- 
lytical approaches developed in |jlj, [l?], |o| . It should 
be stressed also that our results have relevance to a wider 
class of quantum systems, like e.g. quantum spin glass or 
other disordered systems with interactions. Indeed, the 
existence of highly degenerate classically stable configu- 
rations is a general property of such systems. 

The paper is organized as follows. In Section || we out- 
line the model, its quantum features and the main points 
of our numerical approach. In Section III we study el- 
ementary excitations of the chain at different values of 
h and show that there exists a structural rearrangement 
of the excitation spectrum at certain h—h c . A detailed 



analysis of this rearrangement, given in Section [TVj, in- 
dicates that we have a quantum phase transition. Our 
results are summarized in Section [V]. 



II. THE QUANTUM FRENKEL-KONTOROVA 
MODEL 

A. Definitions and outline of quantum features 

The model describes a one-dimensional chain of parti- 
cles with harmonic couplings placed in a periodic poten- 
tial. The Hamiltonian reads: 
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p? 

2m 



(3 cos(xi/d) 



a) 



where s is the number of particles, Pi and their 
momenta and coordinates, and in the quantum case Pi = 
—ihd/dxi. In this paper we use units m = d = a = 1, 
that correspond to dimensionless Hamiltonian (see for 
details e.g. 
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K cos(xi) 



(2) 



with a dimensionless parameter K = (3 /ad 2 and the di- 
mensionless Planck constant H is measured in units of 
d 2 \fmoi. In this way, K is the chaos parameter in the 
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Chirikov standard map |23| , We use the standard 
boundary conditions 



xq = 0, x s = L . 



(3) 



where the chain length L = 2-7T ■ r consists of r peri- 
ods/wells of the external field. We analyze the stan- 
dard case of golden mean ratio corresponding to r/s — > 
(-y/5- l)/2 (see, [TJ, The potential energy of 

the chain 



U({x}) = J2 



l) 2 



if cos(a:i) 



(4) 



has a large number of minima, corresponding to differ- 
ent possible distributions of particles among the wells. 
The classical ground state (absolute minimum of the po- 
tential energy) is characterized by some special ordering 
discovered by Aubry jf], |23|, |24j. In addition, there are 
also local minima, known as "configurational excitation" 
states|Q. Many of them have exponentially small en- 
ergy separations from the energy of the classical ground 
state In the classical case all these states are well 
defined (distinguishable) and absolutely stable. 

In contrast, in the quantum world these states are 
metastable due to quantum tunneling. Since they are 
practically degenerate with the classical ground state the 
actual quantum ground state is built by many of them. 
The admixture of metastable classical configurations in 
the quantum ground state was first discussed in jfj], [FJ|. 
Let us summarize below the most important aspects of 
the quantum ground state, related to the tunneling be- 
tween these configurations. 

In the quasiclassical region, the transition amplitudes 
between different metastable configurations are exponen- 
tially small. Hence, the ground state wave function is 
(with exponential accuracy) a sum of non-overlapping 
parts, each referring to a particular classical configura- 
tion. Any average over the quantum ground state is 
(within the same accuracy) a weighted sum of averages 
over relevant classical configurations. Note that in this 
limit the main contribution to quantum motion comes 
from phonons which characterize small vibrations around 
a classical equilibrium configuration. In fact, the phonon 
spectrum is only weakly sensitive to the choice of specific 
configuration |3lJ and the influence of tunneling processes 
on the global (thermodynamical) properties of the chain 
is expected to be small. In particular, the phonon spec- 
trum is still characterized by a phonon gap similar to the 
classical case. In this regime the tunneling transitions be- 
tween metastable configurations are very slow comparing 
to phonon frequencies and can be considered as well sep- 
arated instantons. 

At higher values of h the tunneling rate increases 
and becomes comparable with the frequency of phonons. 
In this regime phonon oscillations are large and essen- 
tially anharmonic and therefore the known analytical ap- 
proaches Jl4| , [H| [Tj| Q based on the gauss-like wave 



function profile are not quite adequate. Instead, inter- 
actions between instantons and phonons come into play 
here. As it will be shown later, this leads to a new sliding 
phase appearing at h > h c . 

In fact, the first effects of quantum tunneling between 
classical configurations were already seen in the numeri- 
cal studies [jrj[ (see Fig. 4b). However, they were not at- 
tentively analyzed and the mixture of classical metastable 
configurations induced by quantum tunneling was not 
discussed. Certainly, tunneling processes affect signif- 
icantly the average positions of particles at large time 
scales. Therefore the quantities based on mean expecta- 
tion values of particle positions are not adequate even in 
the deep quasiclassical regime. Hence the basic concepts 
of classical treatment, like hull J|, |23|] and g- functions 
|PH , [l2|| , are also not adequate in this quantum regime. 
The point is that in the case of particle tunneling be- 
tween two classical equilibrium positions the expectation 
value does not correspond to any probable (in the clas- 
sical sense) particle position. This is evident from the 
analogy with the two-well potential problem, where the 
mean expectation value of a particle position coincides 
with the top of the barrier. 



B. Path integral and numerical simulations 

Tunneling effects are best understood in the Feynman 
path integral formulation of quantum mechanics |Q. In 
particular, the transformation to "Euclidean" time vari- 
able r: t — it makes the tunneling transition, or in- 
stanton, local in time r. In the quasiclassical regime 
the tunneling probability is very small and therefore the 
mean separation between instantons is large compared to 
their size. This allows to use the approximation of dilute 
instanton gas |7J]. 

Feynman path formulation in the "Eucleudean" time 
allows direct numerical simulations of quantum systems 
with probabilistic treatment J33f of the path integral: 



1 



Dx[T]exp{--S[x(T)}) 



(5) 



where the "Eucledian" action 



S[x(t)] 



4^ V 2 



K cos{xi) 



(6) 

has the same form as a total energy of a chain of particles 
in the usual real time variable, integrated over some time 
interval [0, To]. This is equivalent to the consideration of 
the quantum system at some finite temperaturep3j 



T = h/r 



(7) 



We assume periodic boundary conditions in the time 
direction, which correspond to a path closed on a torus 



Xi(ro) =ii(0), £ = l,...,s - 1, 



(8) 
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and integrate over the initial conditions Xi (0) in order to 
restore homogeneity of paths along the time torus. This 
allows to improve the statistics in data measurements by 
averaging the data along the torus. 

The numerical simulation of the path integral (||) 
needs discretization of the time variable t„ = At • n, 
n = 1,...,N by splitting the time interval To into N 
steps of size At = tq/N. As a result the original time- 
continuous model (||) turns into a 2-dimensional lattice 
model with the action 



N s 



s = EE 



2Ar 



+ 



+At At • K cos x^ n \ . (9) 

The time step At should be chosen small enough to 
approach the continuous limit of the original model (||) . 
This leads to the following requirements. At any time 
step At the p ath variable Xi jumps by a random shift 
Axi ~ VhAr. The first obvious requirement is that the 
shift Axi should be small compared to the spatial scale 
of the potential energy variation. Another requirement 
comes from the standard derivation of the path integral 
| ^3j |34]| : the potential energy terms should be small com- 
pared to the kinetic energy: 



At- 



{x. 



i+l,n 



,r 2 



K cos Xi 



< 



n+l 



X- S 1 



2Ai 



(10) 

The latter condition makes our lattice anisotropic with 
respect to spatial/time directions. In spite of the com- 
mon belief [35| on the equivalence between quantum 1- 
dimensional chains and classical 2-dimensional statistical 
(e.g. XY) models, the above discussion shows that some 
particular care has to be taken. 

An accurate treatment of particle tunneling between 
two wells in the chain requires some modifications of the 
action (Q), For a finite (not very small!) time step At, 
there exists a probability for a particle to jump over the 
potential barrier in one time step|Q. The corresponding 
path contribution to the action will be strongly underes- 
timated by eq.@ and, respectively, the probability of 
tunneling will be too high. To cure this problem we use 
in our simulation an improved |34|] version of action (||) , 
with the potential energy term AT-U(x n ) replaced by the 
integral At • J x " +1 U(x)dx/(x n +i ~ x n ) along a straight 
line that links subsequent (in time) points {x n ,x n +i) of 
the particle path. This significantly improves the accu- 
racy of numerical simulations. 

For simulations of path ensembles we use the standard 
Metropolis algorithm j3(J. Each iteration looks as fol- 
lows: at any fixed time slice at number n we update 
sequentially the particles coordinates a^n, i = 1, ■ ■ • , s; 
then we go to the next time slice: n — ► n + 1, and so on. 

The system has two principal characteristic relaxation 
time scales that are originated by different underlying 
processes. Basically one can estimate the number of iter- 
ations as Nn ~ l/uj( min >AT, where o;( mm ) is the lowest 



frequency relevant to the process and At is time dis- 
cretization step. The shortest scale is related to the path 
relaxation with respect to main phonon modes, and is 
typically of the order of a few tens of iterations in the 
quasiclassical regime, where phonons have a gap of or- 
der unity. Another time scale is much larger and is de- 
termined by the smallest frequency related to either the 
tunneling rate between different classical configurations, 
or to the lowest phonon frequency available in the system 
for h > h c . 

For the initial state we choose the Aubry classical 
ground state. Then, applying iterations, we generate a 
path ensemble. To be sure that the system does actually 
relax to statistical equilibrium with respect to the slowest 
processes described above, we control the mean number 
of particle path crosses a top of the potential barrier and 
we discard all configurations in the ensemble until this 
quantity stabilizes. All computations are done for the 
chaos parameter K — 5. The required number of itera- 
tions to reach relaxation is very sensitive to the value of 
quantum parameter h: for example at h — 3 this number 
is N lt ~ (2 4)10 2 , while at K = 1 it is of the order 
of 10 5 ^ 6 . This explains why in the first studies @ |l|| 
done at N it < 10 4 many details at ft < 2 were not seen. 
We study chains with up to 233 particles. 



III. ELEMENTARY EXCITATIONS 

The most important information about the quantum 
system is contained in its spectrum of low-lying elemen- 
tary excitations. Being the net manifestation of sys- 
tem internal structure it reflects any structural transi- 
tion which can occur in the system, and it provides a 
complete description of low-temperature thermodynamic 
and kinetic properties. We extract this spectrum using 
a novel approach based on the analysis of Fourier spec- 
trum of Feynman paths. This approach provides a most 
direct way to see and resolve different excitations in the 
system. Then we compare our method with the more 
traditional one, based on the study of time correlation 
functions. This comparison provides a self-consistency 
check and demonstrates the advantages of our method. 



A. Spectral properties of Feynman paths 

In classical nonlinear dynamics, the Fourier analysis of 
trajectories plays a key role in understanding of periodic 
motion of complex systems. In a similar way, the spectral 
characteristics of Feynman paths are closely related to 
the properties of elementary excitations in quantum sys- 
tems. We start our studies from the quasiclassical limit 
h —>■ 0, where this relation is exact, and extend them to 
higher values of h. 

Let us consider the Fourier image of the path variable 
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Xi(r): 

1 f T0 

cti(uj m ) = —— / dTXi(r) exp(iu) m T) (11) 
\To Jo 

where u> m = mzu, w = 2it/tq, and — oo < m < oo. 
The path variable Xi(r) is real, therefore a.j(— u> m ) = 
(oj(a; TO ))*. To get insight into the physical content of this 
quantity let us consider the quasiclassical regime 
Then, for small variations Xi(r) = X{ + foj(r) around 
the classical static trajectory {x{\, one can expand the 
action (g) up to second order terms in 8xi (r). Next, us- 
ing the spectral expansion for <5iEi (r) and performing the 
integration over t, one gets: 

S = S [x} + 

+ ^2^2 2^ LJ " lSlk ~ n fk)a-i{~uj m )a k {uj m ) (12) 

m i,/c 

^fc = *ifc(2 + if cos(aii)) — <y*,fe_i -5»-i,fc (13) 

Now, by the transformation to normal modes A^(w m ) = 
J2i ^''ii( w m)i where eigenvectors V^'' 1 satisfy the equa- 
tion fiffcl^ = vfV^ one gets the standard representa- 
tion of the action as a sum of independent phonon modes: 



m i 



(14) 



Finally, the path integral (||) turns into a product of or- 
dinary integrals 

Z = JJdAW(0) exp — z/ ; U W (0) /2ft) x 
3 i=i ^ ' 

oo 

Y[ dReA {l) (uj m )d\mA (l) (uj m ) x 

m— 1 

and one arrives to the well known result for the correlator 
of free phonon modes: 



(A®(«j m )A«>(w m ,)) = 



(15) 



We note that this result is obtained from Eq.(|6|) where 
the action is continuous in time variable. In the dis- 
cretized version (JoJ) the number of harmonics is fi- 
nite: \m\ = 0,1,..., M, where M = t /2At. It can 
be shown that the only modification induced by dis- 
cretization is the replacement in (|l5|): uj m — > oj m = 
(2u>m /it) sin(7ro; m /2a;M), % = Mm. As a result, the 
spectral function for phonons is given by 



F«(5„0 = < A«(w ro ) 



(16) 



Hereafter, instead of LJ m we assume its discretized version 
uj m , and in the following the tilde will be omitted. 

The expression ([l5]) is the well known Wick rotated 
Green function (in the frequency representation) for a 
single free particle in the phonon field theory, which has 
in our case one spatial dimension. 

In the quasiclassical regime the amplitudes of phonon 
oscillations are small and the interactions between 
phonons due to anharmonicity of the Hamiltonian (||) 
are negligible. At higher h, the amplitudes of phonon 
vibrations grow as h 1 / 2 and their interactions become 
more important. In general, interactions can essentially 
modify the Green function ( fl5| ) for phonon excitations. 
This actually happens for h > h c where the spectrum of 
excitations is significantly changed. However, one may 
expect that the spectral function of elementary excita- 
tion remains of the same form 



(17) 



which differs from (|T^) by the renormalized frequency 
value vi and by an overall renormalization factor /, in 
analogy with the Green function behavior in renormaliz- 
able quantum field theories (see, e.g. j29|). In fact, this 
idea is well supported by numerical data. 

An extended elementary excitation involves all parti- 
cles in the chain. In turn, the Fourier harmonics of any 
particle coordinate in the chain are a sum of contribu- 
tions of many elementary excitations 



ai(u) m )\' 



0) 



(18) 



where the sum goes over all chain excitations. In the deep 
quasiclassical case the main contribution comes from 
phonons modes, and the sum goes over phonon modes 
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V; 



(I) 



and 



vi- 



The goal of our study here is a complete picture of low- 
lying quantum excitations in the chain, both for low and 
high values of h. In fact, the spectrum of low-lying excita- 
tions is crucially dependent on H: there are domains with 
a qualitatively different behavior. In order to illustrate 
this, let us consider the amplitude of quantum motion 
of particles in the chain, given by Eq.([l8|). It is seen 
that the contributions of low-frequency modes (small v{) 
dominate in the limit uj m — > 0, provided that their wave- 
function profiles are not small at the particle position 
i. Therefore the spectral function (|lgl), computed in this 
limit, gives a rough estimate for the frequency v\ of the 
lowest mode. 

In Figjl] the amplitude of the lowest Fourier harmonic 
aj(oJi) with u>i = zu = 2it/to is plotted as a function of 
the particle position i at different h = 0.6 — 8. One can 
see that the whole interval of h splits naturally in three 
regions of qualitatively different behavior: (i) the qua- 
siclassical region K < 1 where the amplitudes a^c^) of 
the harmonics are very small and depend on the particle 
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In I d: 




harmonic form 



Figure 1: Dependence of the amplitude of the lowest Fourier 
harmonic 04 = aj(ui) on the particle position i in the chain 
at different h. Here uj\ = zu = 2ty/to, m = 1. The chain 
parameters are s/r — 89/55, K = 5, to = 80. Typical number 
of iterations is (1.5 -7- 5) • 10 5 at each value of K. 



positions in some regular way; (ii) the transition region 
1 < fi. < 2, where this dependence is highly irregular and 
interactions between instantons and phonons are impor- 
tant; and (iii) the region h > 2, where this dependence 
becomes regular again and where, as we shall see be- 
low, a new phonon branch appears. Let us note, that at 
H < 1 the regular structure along the chain is quasiperi- 
odical, which reflect a fact that a classical chain is built 
of "bricks" of two principal sizes |28| . Above Sw2 bricks 
are "melted" and chain properties become even more ho- 
mogeneous along the chain. In the intermediate region 
irregular peaks come from different non-overlapping in- 
stantons contributions, which as any tunneling effects are 
highly sensitive to small variations of potential barriers. 
Below their contributions are exponentially small, and 
above, as we see further, they overlap and form new slid- 
ing phase of the system. 

In the following we analyze in detail these regions, cor- 
responding to different intervals of ti. 



B. 



Quasiclassical region h < 1 



For H < 1 the tunneling between different metastable 
classical configurations is negligible, and the particles 
mainly vibrate around some classical equilibrium po- 
sitions. In this case the elementary excitations are 
phonons, and the quantization of the chain is reduced 
to the quantization of phonon modes, see e.g. Jl8| . 

To single out low energy excitations we use the follow- 
ing approach. Of course, a particular excitation can be 
selected if the corresponding mode Vj 1 ^ is known, but in 
general this is not a trivial task. However, for low lying 
excitations one may expect that the modes have a simple 



V s {k t ) = V, 



(') 



sin(fcij), 



ki = ttI/L, (I = 1,2,...) 



(19) 



where the wavelength Aj = 27r/fc; is much larger than 
a characteristic size of inhomogeneity in the chain. Di- 
rect numerical computations |31[ of phonon modes in the 
classical FK chain support this ansatz (fTs|). 




Figure 2: The phonon spectral function F^(ui m ) versus the 
rescaled frequency Lu m , data are shown for the lowest spatial 
mode with / = 1. The chain parameters are: s/r = 89/55, 
K = 5, ro = 80, h — 0.8. Here m varies from 1 to 89, but for 
clarity only selected values are shown. The solid curve gives 
the fit by Eqs.( |l6|]l7| ), open circles show numerical data. The 
fit determines the phonon frequency of the first mode (v'{ = 
3.170 ± 0.012) and the renormalization factor / = 1.0034 ± 
0.0034. 

The numerical test of this anzatz (|l^) is given in 
Fig.||. Here a typical result of quantum simulations of 
the spectral function F^(u m ) is shown for the lowest 
phonon mode I = 1 at K = 0.8. Fitting the data by 
Eq.(p"7|) with the renormalization factor / and the fre- 
quency of the phonon mode vi as free parameters, we 
obtain / = 1.0034 ±0.0034 and v\ = 3. 170 ±0.012. This 
fit shows that the renormalization factor / is remark- 
ably close to unity, in spite of the fact that the value 
h = 0.8 is not small. Hence, the ansatz ( |l9[ ) provides a 
good approximation to actual profiles of lowest phonon 
modes. This also indica tes th at the Gauss approxima- 
tion used in the Section [II A and in papers (if], |l^, |2c|l 
works fine here. However we note that at the same time 
the quantum effects renormalize substantially the phonon 
frequency (y\ = 3.170) compared to its value at h = 0.2 
{y\ = 3.706 ± 0.019). The computations for the classical 
FK chain give v\ = 3.717 (h = 0). 

These data show that at small h the frequency v\ ob- 
tained from the quantum simulations approaches to the 
frequency of phonon mode in the classical chain. This 
fact gives a very important check of the consistency of 



7 



our quantum simulations. We stress that a good agree- 
ment between the numerical data of FigJ| and the the- 
oretical spectral function ( |l7| ) takes place in the whole 
frequency range of uj m . This is, in fact, a very important 
consistency check of the good relaxation of our paths en- 
semble at all frequencies, including paths fluctuations at 
the lowest frequency available in the system. 
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Figure 3: The phonon dispersion law u{k): open circles show 
the data obtained from the fit as it is shown in Fig.^ for I = 
30, and the straight line shows the fit given by Eq.jiol). The 
chain parameters are the same as in Fig]!], the wave vector 
k = nl/L. 

By fitting the data for different phonon modes I = 
1-^30 one can extract the dispersion relation for phonons 
v(k) where k = nl/L (see FigJ|). The majority of data 
points follow the straight line given by the formula 



u 2 (k) 



2 2 ; 2 

i/n + c k 



(20) 



where vq is the phonon frequency gap, and c is the 
velocity of sound. The fit of numerical data gives: 



= 3.204 ± 0.026, c 
3.706 ±0.019, c 2 = 



= 15.0 ± 0.7 for h = 0.8; 
13.3 ± 0.5 for h = 0.2. These 
quantum data should be compared with the classical case 
where v 2 = 3.697, c 2 = 11.5 (H = 0). We note that there 
is a difference between the frequency of the first spatial 
harmonic vi and the frequency gap value vq obtained 
from the dispersion law. However this difference is small 
and comparable with the statistical errors. For small h 
the parameters of the dispersion law converge to their 
classical values. 

The described approach allows to obtain a complete 
information about low-energy phonon excitations in the 
whole quasiclassical region h < 1. 



C. Transition region 1 < h < 2. 

As it is seen from Fig.0, this region corresponds to 
the transition between two regimes h < 1 and h > 2 



where the dependence of the quantum excitations on the 
particle location in the chain looks quite regular. 

As it will be shown later the irregular behavior in the 
region 1 < h < 2 is related to a significant increase of 
the density of instantons. At high density the interac- 
tion between instantons becomes important and results 
in onset of new phonon branch at h > 2. In this section 
we discuss the properties of instantons and phonons and 
obtain estimates for their frequencies. 

Let us start with a discussion of tunneling effects. 
For h < 1 the contribution of tunneling to the spec- 
tral function is exponentially small being proportional 
to exp(— const /K). However at H > 1 the tunneling prob- 
ability becomes large and it gives a significant contri- 
bution to the spectral function. The transition to this 
regime is seen in Fig.|j] as a sequence of sharp isolated 
peaks. Following the pioneering paper J32[ , a tunneling 
event can be associated to an instanton. In the imagi- 
nary time representation, the instanton is a local jump 
between two wells, which is fast compared to the mean 
time interval between subsequent jumps: while the size 
of instantons (in time r) is practically independent of h, 
the separation between them is exponentially large in the 
quasiclassical limit (low instanton density). In particu- 
lar, this means that in the first approximation one may 
consider instantons as independent jumps, as can be also 
checked from a direct examination of our path ensem- 
ble. Here we should stress two important properties of 
instantons in the quantum FK chain: 

(i) each jump of a particle i in its position Xi gives 
displacements of neighboring particles, which decay ex- 
ponentially with the distance from the jump location, i.e. 
instantons are exponentially localized in space inside the 
chain; 

(ii) instantons are distributed inhomogeneously along 
the chain since the tunneling probability is highly sensi- 
tive to variations of barrier heights due to chain inhomo- 
geneity. 

A simple explanation of the exponential localization of 
an instanton (along the chain) comes from the fact that 
the new static configuration produced by it can be seen as 
a local static defect on the original configuration, which 
is known to dye away exponentially with the classical 
Lyapunov exponent || |2^, Indeed, in Fig{j] one can 
see that at H ~ 1.2 -j- 1.3 instanton contributions are 
exponentially peaked around some particular positions 
along the chain. 

Let us now consider the properties of elementary exci- 
tations originated by instantons. There is a question how 
to select numerically a single instanton excitation. Obvi- 
ously, the ansatz C P- 91 ) used for phonons is good only for 
extended modes while instantons are localized in space. 
Therefore in this case we analyze numerically the fre- 
quency spectrum of a given particle i in the chain (de- 
fined by Eq. (|ll])) . If h is not too high (close to one), 
instantons do not overlap, and the main contribution to 
the spectrum comes from the instanton which is near to 
the given particle. This contribution reaches its maxi- 
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mum for a particle which actually jumps. 

At the same time besides instanton jumps the quan- 
tum motion of a particle in the chain contains a contribu- 
tion of many phonons with different frequencies (see, e.g. 
Eq.(|l8|)). This phonon background should be subtracted 
in order to single out the contribution of instanton. For- 
tunately the frequencies m of phonon excitations are 
much higher than the frequency zA mst ) of chosen instan- 
ton. Hence these two types of excitations are well sepa- 
rated in the frequency domain. Therefore in our analy- 
sis of the instanton contribution to the spectral function 
Fi(uj m ) = (|<Xj(u> m )| ) we can restrict ourselves to the fre- 
quency domain uj m < tutound, where the boundary unbound 
is chosen by the condition j/ mst ) <§; uJbound <C v\ for all 
phonon modes I. In this frequency domain the phonon 
contribution ( |l6| ) is practically independent of u> m and 
can be replaced by some constant Cj. Thus we can ex- 
tract the frequency y( mst ) and the weight fy nst ^ f the 
instanton excitation by following fit for the spectral func- 
tion Fi(u> m ); 



fiW = ft st) N fa + (vf nst) Y) + a. (2i) 



(inst) 



ff nst) and 



This fit contains three free parameters v, 
Ci, 

If instantons do not overlap then one may expect that 
there are groups of particles which motion is dominated 
by a single instanton. Inside each group the fit (|2 

should give the same values for the frequencies v\ 

while the variation of weight f^ m with i determines the 
instanton profile along the chain. This case is illustrated 
in Fig.||a, which corresponds to the early onset of instan- 
ton contribution at h = 1.2. The six peaks in the bot- 
tom part of Fig.^ja show six non overlapping instantons, 
while the top part shows the corresponding frequencies 
as a function of particle index i inside the chain. The 
peaks have different amplitudes, and the highest three 
of them involve groups of three particles which have the 
same frequency inside each group. 

At higher H the number of instantons starts to grow 
rapidly and they begin to overlap. A direct confirma- 
tion of this trend is seen in Fig.|]b which corresponds to 
h = 1.8. Here all instantons have about ten percent over- 
lap with their neighbors and their interaction is rather 
strong. As a result the step-like structure of frequencies, 
seen at the top of the Fig.^a, is practically destroyed. 
Thus the instantons are "collectivized" and the consid- 
eration of a single instanton as one-particle jump over 
barrier becomes not adequate. At higher h values this 
process leads to appearance of a new phonon mode. 

Let us now discuss the phonon properties i n the re- 
gion 1 < h < 2. As in the previous Section MB , we 
extract them from the the spectral function F(uj m , k{) = 
F®((J m ) obtained on the basis of the ansatz ([[9]) for a 
phonon mode I. However, in this region of h the spec- 
tral function F^(uj m ) has an admixture of instantons, 
which grows rapidly with h: a change of h from 1.1 
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Figure 4: Dependence of the instanton frequency v, 
and its weight J^" 13 ') (bottom) on the position i inside the 
chain (see text for explanations), (a) h = 1.2. Instantons 
do not overlap: the three highest peaks involves a group of 
three particles each, (b) h = 1.8 Instantons overlap. Chain 
parameters are s/r = 34/21, K = 5, To = 320. 
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Figure 5: Rescaled spectral function F R (uj,ki) = 
(yx/K)F{u), ki) at the very beginning of instanton onset. The 
bump which appears at u < 0.2 corresponds to instanton ad- 
mixture to the phonon spectral function. Black points, open 
circles and stars correspond to h = 1, 1.1 and 1.2. Lines show 
the fit to numerical data with Eq.(J22j). Chain parameters 
used in simulations are: s/r = 34/21, K = 5, to = 320. 



to 1.2 results in more than ten times of the admixture 
weight(see Fig.||). We note that in the absence of in- 
stanton contribution the rescaled phonon spectral func- 
tion F R (uj,ki) ee (v?/h)F(w,h) plotted in Fig.| should 
have an universal limit equal to unity independent of the 
value of h. Hence the increase of F R (ui,ki) at small u> 
stresses the important contribution of instantons. These 
instantons have different frequencies and their contribu- 
tion to the spectral function can be rather complicated. 
However, we can use again the strong frequency separa- 
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tion between instanon and phonon excitations. Indeed, 
for uj m > i/j » v ( inst ) ) a ll instanton contributions have 
an universal behavior oc w~ 2 . Therefore we may replace 
them by a single "instanton contribution" with some aver- 
age instanton frequency y mst . Then the spectral function 
can be fitted by a sum of two contributions: 

F(u m ,k t ) = f ph h/ (iul + vfih)) 

+f inat h/ (u^ + (9 inst (h) 2 )) (22) 

where f P h(k), vi{k) and fi n8 t(k), v m (k) are free fit param- 
eters for phonons and instantons, respectively. 

Fitting the data for different phonon modes at I = 
1 — 30 we extract the phonon dispersion law (|2(i| ) (see 
Fig.|6|, compare with fit procedure for h < 1). Contrary 
to the case of Fig.||at h < 1, now the data for the disper- 
sion law u(k) are scattered inside some finite band. This 
indicates that the anzatz ( |l9| ) for the phonon profile is 
not so good to single out particular phonon modes. Ac- 
tually, the width of the band provides some measure of 
the inaccuracy. Nevertheless, the phonon modes are still 
approximately defined in the domain of H under consid- 
eration. Their frequency decreases as h — > 2 but remains 
separated from zero by a finite gap. 



4.5 




Figure 6: Frequency vW) of phonons versus wavenumber k 
obtained from the fit (E2) at h = 1.5 (in the middle of the 
transition region). Chain parameters used in simulations are: 
s/r = 89/55, K = 5,r = 80. 



We also note that the fit (|22|) allows formally to deter- 
mine the dispersion law Vi ns t{k) for the instanton branch. 
However the numerical data give irregular scattering of 
points inside a band < vf nst < 0.01 without any clear 
dependence on k. The reason of such behavior is sim- 
ple: projections of irregular positions of instantons on 
the harmonic ansatz ( |1S| ) produce random weights for 
contributions of different instantons. This result repre- 
sents another manifestation of glass-like structure formed 
by instantons frozen/pinned inside the chain. Since the 
positions of instantons are random the phonons cannot 



propagate along the chain on large distances. In fact 
they become localized by disorder in a way similar to the 
one-dimensional Anderson localization (more details on 
the phonon properties in this regime will be presented 
elsewhere ||3l|| ). 



D. 



New sliding phonon branch at ft > 2 
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Figure 7: A sample of quantum paths of particles inside some 
chain fragment, which corresponds to periods of the external 
potential with numbers 22 ~ 30. Dashed lines show bottoms 
of the wells, thick solid lines show the tops of the barriers. 
Note an example of highly correlated instanton transitions 
at r = 15 — 30 which involve particles in up to 5 periods of 
the potential. Parameters of the simulation are: h = 2.3, 
s/r = 89/55, K = 5 and to = 80. 

From Fig.|| one can see that the variation of the ampli- 
tude ai(u>{) with i (low frequency excitations) becomes 
rather smooth at h > 2. This means that the instantons 
are strongly overlapped here. In fact, this means that 
several particles of some chain fragment jumps from one 
well to the next simultaneously (see Fig.0) . This is noth- 
ing but sliding of a local chain fragment along the period- 
ical potential. The typical size of such fragments should 
grow with h. If their sizes reach the size of the chain then 
the sliding mode becomes open, and the phonon gap dis- 
appears. At this point the pinned instanton glass turns 
into the sliding phonon gas. 

A confirmation of this picture is presented in FigJ|, 
which corresponds to h = 2.5 being just above the tran- 
sition point h c w 2. The numerical data for the disper- 
sion law u(k) in Fig.^j are obtained from the fit (22) at 
h > 2. Here the behavior of phonon and instanton modes 
changes dramatically: phonons data v(k) are now irreg- 
ularly scattered over a wide band while the data points 
for instanton branch follow a single line, reproducing 
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Figure 8: The dispersion law v(k) for sliding phonons at 
h = 2.5. The numerical data (circles with error bars) are 
obtained from the fit ( [i^ ) for the instanton branch u inat (ki). 
The straight line shows the best fit ( po| ) to numerical data. 
The chain parameters used in simulations are: s/r = 89/55, 
K = 5,to = 80. 



fairly well a phonon-like dispersion law with zero gap (see 
Figs.|,§). In particular, the fit @ gives z/ = 0.04±0.01, 
which is close to zero. In fact, this value is smaller than 
the minimal frequency 2tt/tq (s» 0.079, at To = 80) and 
therefore it is compatible with zero. 

On the contrary, the wide scattering of data points 
for the phonon branch indicates that the ansatz ( |l9| ) is 
not good for phonon contribution at ft > 2 (see Fig.^J) . 
This scattering of points is related to the localization 
of high-frequency phonon modes. In contrast, a smooth 
behavior of data points for the instanton branch demon- 
strates that the instanton wave functions become close to 
the harmonic wave ansatz (|l9|). Hence these excitations 
are delocalized. This leads us to the conclusion, that for 
ft > 2 the instanton branch is replaced by a new gapless 
branch of new sliding phonons. 



E. Global picture of elementary excitations in the 
FK chain 



The ensemble of data for the dispersion law v{k) of 
elementary excitations at different values of ft is shown in 
Fig.^ by the two sheets representing the phonon branch 
(1) and the instanton branch (2). The numerical data 
are obtained with the ansatz (|l|) by the fit (^). 

The sheet (1) refers to phonons originated from classi- 
cal phonon modes which are well reproduced in the limit 
ft — > 0. The frequencies of these modes are well sepa- 
rated from zero by a large gap. Therefore one may say 
that they form the optical phonon branch. At ft > 1 
these modes show a tendency to become softer, and at 
ft > 1.5 their dependence on k and ft becomes irregular. 
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Figure 9: The frequency v(k) of elementary excitations versus 
the wavenumber k at different h. Chain parameters used in 
simulations are: s/r = 89/55, K = 5, to = 80. The two 
sheets (1) and (2) refer to phonon and instanton excitations 
respectively. 



As it was explained in the previous section this irreg- 
ular dependence is related to the glass-like structure of 
randomly pinned instantons which density increases with 
the growth of ft. The sheet (2) appears at ft < 2 from 
the instanton contributions into the Feynman path inte- 
gral. For ft > ft c m 2 this instanton branch turns into a 
new gapless branch of sliding phonons (see the discussion 
related to FigJ|). 

A more quantitative picture can be obtained from the 
numerical data for the gap i>q and the sound velocity c 
for both sheets in Fig.||. For the optical phonon branch 
the values of v$ and c are determined from the fit ( |2^ ) 
for different values of ft (see Fig .[!(]). The data show that 
the phonon gap remains finite and large at ft < ft c » 
2. In contrast, the sound velocity c drops to zero as ft 
approaches the value ft c . This decay is compatible with 
the fit 



a(ft c - Kf 



(23) 



shown by the solid line with a = 20. 5± 1.3, ft c = 2.0±0.1 
and the critical exponent a = 0.52 ± 0.07. Even if the 
numerical data for c have certain fluctuations they still 
clearly indicate the quantum phase transition at ft c w 
2, where the sound velocity c drops to zero. Further 
extensive numerical studies are required to determine the 
behavior in the vicinity of transition in a more precise 
way. 

For the sliding phonon branch we extract the parame- 
ters vq and c from the data of sheet (2) in Fig|9 using a 
more general fit given by 



v\k) 



+ c 2 k 2 /(l + k 2 /k 2 B ). 



Compared to the standard case (20) we introduce an ad- 
ditional parameter k B to take into account the saturation 
of v(k) at large k (see, Fig.|||}]). The numerical data show 
that the gap v$ is small and does not exceed the minimal 




Figure 10: Dependence of the phonon gap vo and sound ve- 
locity c on ft for the case of Fig.gl at K = 5. Open symbols 
correspond to the sheet (1) in FigM Circles and triangles are 
obtained at to = 80 for chain sizes s/r = 89/55 and 233/144 
respectively. Stars correspond to s/r = 34/21 and to = 320. 
The solid line (bottom) gives the fit for 1.3 < ft < 2 
(see text). The full circles (top) for ft > 2 refer to the sliding 
phonon branch from the sheet(2) of Fig.|i], they indicate a zero 
phonon gap. 



frequency in the system zu = 2it/tq. Hence, the gap vq 
is compatible with zero (see Fig|l(] (top)). At the same 
time the sound velocity c for the sliding phonon branch 
grows approximately linearly with h (see Fig Jill). The 
best fit of numerical data for ft > 2 gives 



a ■ (h — b), 



with a = 7.6±0.1, b = 1.57±0.05. Formally, the value of 
b is different from the value H c — 2.0 in (p3|). However in 
view of large statistical fluctuations both fits for optical 
phonons and sliding phonons are compatible with the 
quantum phase transition at h c w 2. 



F. Time correlations 

Time correlations are closely related to the frequency 
spectrum of elementary excitations in the system. Their 
analysis is probably the most traditional way to ex- 
tract properties of the elementary excitations, see, e.g. 
|?3L M . Below we discuss the connection of this tradi- 
tional method with our approach based on the Fourier 
spectrum of Feynman paths described above. We com- 
pare the results obtained by these two different methods. 

The one-particle time correlators are directly related to 
the Fourier spectrum of Feynman paths, and are given by 



Figure 11: Dependence of sound velocity c on ft for excitations 
on the soft phonon sheet (2) in Fig.[j. Points correspond to 
chain parameters s/r = 89/55 K — 5 and to = 80. The solid 
line shows the linear fit to the sound velocity data inside the 
region ft > 2 (see text). 



the following expression: 

(xi(O)xi(r)) = — != V" |aj(w m )| 2 exp(- 

. / Tr, ' ■ 



In fact, many elementary excitations with different fre- 
quencies contribute to a particle motion in the chain. In 
order to single out a particular phonon mode we study the 
correlators of normal modes X 1 (t) = ^ V^Xi(r) where 
V® are eigenvectors defined in Eqs.(|l3|)-(|l3). Then from 



(24) Eq.© we obtain 



(X 1 (0)X l (r)> = -L V | A 1 (u m ) | 2 exp {-iLo m r) 



(X l (0)) 2 )(e-^ 



+ e 



— fi(no— t) 



(25) 



where (^(X'(O)) J = h/2vi is the contribution of a single 

phonon mode and the periodicity along the time torus 
results in a second exponential term in (25). 

However, Eq.(|l5|) assumes only small quantum fluctu- 
ations (~ h) around some classical trajectory. But due 
to tunneling effects (or instantons) a particle jumps from 
one classical trajectory to another and its actual mo- 
tion is given by a sum of phonon xf h ^ (r) and instanton 
x (*nst) ^y-j con t r ibutions: 

In general, both motions influence each other, but in 
the quasiclassical limit h — > they can be considered 
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as independent. In this limit they have quite differ- 
ent frequency scales: the phonon frequency v\ is of the 
order of K 1 ! 2 while the frequency of tunneling jumps 
Vi ns t is exponentially small. In contrast, while the am- 
plitude of phonon oscillations in the limit H — ► is small: 

/^Ei^f/r)^ \ oc h, the amplitude of jumps is defined 



by difference between equilibrium particle positions in 
two neighbor wells (in our case it is ~ 3 4), i.e. it 
does not depend on h. Hence, for a jumping particle 

! " \ -10 which is not small even in the 



quasiclassical limit % — > 0. Therefore the instanton con- 
tribution to the time correlator has the form |E3 



(26) 



The pre-exponent factor is large for jumping particles 
even in the deep quasiclassical regime. 

In fact, not any particle can easily jump from one well 
to another: different classical trajectories have different 
actions and all jumps which result in a large change of ac- 
tion AS > h are inhibited. In particular, it is clearly seen 
in Fig]l] that the number of instanton peaks is smaller at 
smaller h since the contribution of transitions with large 
difference in action AS is suppressed. 

Finally, the general form for the time correlator which 
takes into account the instantons contribution takes the 
form: 



ee (X l (0)X l (r)) 

= ((x l (o)) 2 

+ E 



Wi. 



(27) 



Here the first term describes the phonon contribution and 
in the second term the sum is taken over instantons and 
the instanton weights w inst = J2i,k V i^ V k^ ( x i(ty x k(0)) 
describe the overlap with the ansatz (|l9|). Both types of 
contributions ( p7j ) are clearly seen in Figjl^ at h = 1.4. 
The initial rapid drop at r < 1.5 corresponds to phonon 
contribution while the slow decay at r > 1.5 corresponds 
to the instanton contribution. This initial drop is related 
to the existence of large quasiclassical gap for phonon 
excitations. For h — 2.5 > h c the gap disappears and the 
correlator decay very slowly (see Fig.|l2l). 

For H > h c we have a new phase where instantons are 
replaced by sliding phonons. Therefore in this regime we 
fit the numerical data for C(r) by Eq.(p7|) with Wi ns t — 0. 
The results for different I allow to obtain numerically the 
dispersion law v(k) shown in Figjl^ (open circles). How- 
ever the accuracy of this data is not so good compared to 
data (full circles) obtained from the analysis of Fourier 
spectrum of Feynman paths described in the previous 
sections. 
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Figure 12: The numerically computed time correlator C(t) = 
(X ! (0) X 1 (t)) for I = 1 at different time separations r for 
H = 1.4 (lower points) and h = 2.5 (upper points). The 
parameters of the chain are s/r — 89/55 K = 5 and tq = 80. 




Figure 13: Frequency of elementary excitations v{k) versus k 
for sliding phonon branch at h = 2.5. Open circles show data 
extracted from the time correlator C(t) for different I (see 
Fig . p^l ), full circles present the results obtained from the fit 
( ^2[ ) for the frequency spectrum of Feynman paths (see Fig.^|). 
The chain parameters are the same as in Fig.|l2|. 



IV. QUANTUM PHASE TRANSITION 

The structual rearrangement of the elementary exci- 
tations spectrum can be related to quantum phase tran- 
sition in the chain from a "pinned" to "sliding" phase. 
However, we would like to stress that in contrast to the 
classical picture p3|, E4) in the quantum case the absence 



13 



of energy gap for excitations is not necessarily related to 
the opening of the sliding phase. Indeed, due to quan- 
tum tunneling through Peierls-Nabbarro barriers related 
to instantons there are excitations with energy which de- 
creases exponentially with the increase of barrier heights. 
Formally, this corresponds to the disappearance of exci- 
tation gap. Therefore, to confirm firmly the appearance 
of sliding phase we need to consider spatial correlations of 
particle motion in the chain. The sliding phase appears 
when the spatial correlation length becomes comparable 
with the length of the chain. 




A. Spatial correlation length 



The analysis of spatial correlations (correlations be- 
tween the motion of different particles in the chain) 
is the most evident way to observe the transition be- 
tween pinned and sliding phases. In principle, the spa- 
tial correlation function can be explicitly computed if the 
spectrum and the wave functions of elementary excita- 
tions are known. In fact, we have a complete quanti- 
tative picture for phonon modes, at least for the low- 
lying ones. To obtain numerically the value of the spa- 
tial correlation length l c we assume that the elementary 
excitation spectrum is given by the dispersion relation 
v 2 (k) = v 2 + c 2 k 2 and the corresponding phonon modes 
are given the ansatz 
correlator reads: 



9h. Then the same-time spatial 



Figure 14: The dependence of the spatial correlation length 
l c on h. Crosses, full and open circles correspond to chains 
with s/r = 34/21, 89/55 and 233/144, respectively; K = 5, 
to = 80. The vertical dot-dashed line marks the quantum 
phase transition at h c » 2. The positions of the horizontal 
dotted lines are proportional to the chain length L — 2nr. 



B. Longwave response 



{(Xi - (Xi)) (Xj - (Xj))) 

cos(k(i-j)) 



Lv(k) 



= —K (\i-j\/l c ), (28) 

7TC 



where Kq(x) is McDonald's function, with a known 
asymptotics: Kq(x) — > x ->oo \J~^ e ~' 1 \ an d lc = c/i/q. 

The fit (|8|) of the numerical data gives the value of l c 
for different values of ft as it is shown in Fig.|lJ. It is seen 
that the length l c has a sharp increase at ft c « 2, and 
for ft > ft c it becomes comparable with the length of the 
chain. This indicates that we have the quantum phase 
transition near ft c w 2. Indeed at ft < 2 the length l c 
is practically independent of the chain length L while at 
ft > 2 it starts to increase with L. This is confirmed by 
the data of Fig. where in spite of strong fluctuations for 
ft > 2 the length l c becomes comparable with the chain 
size L. 



The appearance of the new sliding phonon phase im- 
plies that the response of amplitude of the longwave 
modes should be large in this regime. To test this ex- 
pectation we present the dependence of the amplitude 
A 1 (uj = 0) on ft in Fig.|l5|. The numerical data demon- 
strate a sharp increase of A l (0) near ft c w 2. It is inter- 
esting to note that due to the existence of frequency gap 
for phonon excitations at ft < ft c the amplitude A l (0) is 
not very sensitive to the variations of I. On the contrary, 
for ft > ft c the gap disappears and A 1 (0) starts to depend 
on I (see Fig.|l5|a). In a similar way A l (0) is independent 
of the chain length L for ft < ft c while at ft > ft c it grows 
with L (see Fig.|l5[b). The numerical data of Figjl5|b for 
1.5 < ft < 2 can be described by the fit 



<(^(0)) 2 )«A(ft c -ft)- 



which gives 7 = 5.06 ± 1.72 and ft c = 2.01 ± 0.05. 
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Figure 15: The amplitude square of zero-frequency quantum 
fluctuations (^(O)) 2 as a function of ft for K = 5 and To = 80. 
(a) Full circles, squares and open circles are for modes I = 1, 3 
and 5, respectively, for the chain with s/r = 89/55. (b) Open 
circles, full circles and stars are for s/r = 233/144, 89/55 and 
34/21, respectively; I = 1. The vertical dot-dashed lines mark 
h c 2, the horizontal dotted lines give the average values for 
ft > 2. 



C. Other characteristics 

Another way to test the transition from the pinned 
instanton glass to the sliding phonon phase is to measure 
the sensitivity to small shifts of boundaries of the chain. 
With this aim we consider the shift of boundary particles 
i = and i — s given by 



x (t) = a s cos(2itt/t ), x s (t)=x + L, 



(30) 



where the amplitude as was fixed at as — 0.5. In 
Fig.|l^ we present the dependence of the response func- 
tion R(i) = (xo(xi — {xi))/{xo) on the particle number 
i inside the chain. It is seen that the response in the 
center of the chain drops strongly when the parameter 
h changes from ft = 2.2 to ft = 1.8. This means that 
the chain is locked for ft < ft c 2 while for ft > ft c the 
chain slides following the displacements of the boundary 
particles. 

In order to get a more quantitative picture, we estimate 
the value of the response function R(i) at its minimum in 
the middle of the chain by taking its average value inside 
the central region at i = 12 -j> 22: R m in — (R(i))i (this 
interval is shown in FigJI^ by horizontal dashed lines). 
The dependence of R m in on the parameter ft is presented 
in Fig.[L7|. We see, that the correlator in the central region 
of the chain deviates from zero at ft > 2. According to the 
numerical data the response R m in is very small for ft < 
ft c w 2 while it becomes rather strong for ft > ft c . This 



Figure 16: Dependence of the response function R(i) on the 
particle position i in the chain at ft = 1.8 (full circles) and 
ft = 2.2 (open circles). The horizontal dashed lines show the 
interval of averaging for the minimal response value R m in in 
Fig.^. Parameters of the chain are s/r — 34/21, K = 5 and 
to = 80. 




Figure 17: Dependence of the minimal response R„ 
for the parameters of Fig.|l6|. 



on ft 



confirms the existence of the quantum phase transition 
from the pinned to sliding phase at ft = ft c » 2. 

What is the order of this transition? In our simulations 
we have repeated a cycle changing slowly ft from ft = 1 
to ft = 4 and back about hundred times, but no hystere- 
sis was found within the statistical errors. Indeed, the 
difference of the chain energy per particle at the upward 
and backward paths did not exceed 10~ 3 . This differ- 
ence should be compared with the change of the energy 
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Figure 18: Dependence of the total chain energy per parti- 
cle on h for the chain parameters K = 5, s/r = 89/55 and 
to — 80; eo = 5.302 is an energy point chosen from conve- 
nience/illustration reasons. 



per particle 5E s» 2.4 which takes place when ft changes 
from ft = 1 to ft = 4. The absence of hysteresis excludes 
the phase transition of the first order. We also do not 
see any breaks in the dependence of the chain energy on 
ft. In order to make more visible small deviations from a 
linear law we plot in Fig .|l^ the quantity (e — eo)/K, where 
eo = 5.302. The numerical data show that the slope of 
energy dependence changes near ft = ft c 2. This change 
of slope is located approximately at the same value of ft 
where the divergence of the correlation length takes place 
(see Figjl4|). These data indicate that we have a sec- 
ond order quantum phase transition which appears near 
ft c w 2. However, more extensive numerical simulations 
are required to determine more precisely the order of the 
transition. 



V. CONCLUSIONS 

We have studied quantum tunneling phenomena in 
a particular model of glassy material, the Frenkel- 
Kontorova chain. This system has a lot of states, which 
are exponentially degenerate and (meta) stable in the 
(quasi) classical limit. In the quantum case there are 



tunneling transitions between these states that can be 
understood in terms of instanton dynamics. In the qua- 
siclassical limit at small ft the instanton density is small 
and instantons are local and isolated. With the increase 
of ft their density grows and they start to overlap. 

At sufficiently large ft the instantons are coupled and 
become collectivized. As a result, the tunneling of parti- 
cles in some fragments of the chain proceeds in correlated 
way. The size of these correlations grows until it reaches 
the size of the system. This leads to appearance of ex- 
tended excitations and opening of a new gapless sliding 
phonon branch. Our data show that a quantum phase 
transition takes place between the pinned and sliding 
phases. Absence of hysteresis effects, as well as contin- 
uous dependence of chain energy on ft exclude the first 
order phase transition, so we can classify this transition 
as a continuous quantum phase transition. We stress 
that the quantum phase transition from pinned to slid- 
ing phase takes place in the regime where the classical 
chain always remains in the pinned phase with the finite 
phonon gap. 

The direct analysis of Fourier spectrum of Feynman 
paths ensemble allowed to obtain detailed information 
on the dispersion law of low-lying excitations in both 
quantum phases. Nevertheless, some questions remain 
open for further investigations. For example, one can 
analyze in more detail the effects of interactions between 
instantons at low density and study their propagation 
properties in this regime. Another interesting remark 
concerns the behavior of the system in the vicinity of the 
transition point at ft c w 2. Indeed, in this region the 
kinetic energy per particle is approximately 0.6 that is 
about ten times smaller than the height of the potential 
barrier at K = 5. Therefore more insights are required 
to understand the underlying physics of this transition. 

One can ask on how general are our results obtained 
in the frame of the Frenkel-Kontorova model? In fact, 
we never used any specific features of this model related 
to its nontrivial number theory properties. The only es- 
sential point is the existence of an exponential number of 
quasidegenerate states which is common for glassy ma- 
terials and other disorder systems. Therefore it is very 
interesting to study an analogous quantum phase transi- 
tion in systems with disorder and interactions. 
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